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Although at temperatures T 3> Aqcd the quark-gluon plasma (QGP) is a gas of weakly interacting 
quasiparticles (modulo long-range magnetism), it is strongly interacting in the regime T = (1 — 
3)T'c. As both heavy ion experiments and lattice simulations are now showing, in this region the 
QGP displays rather strong interactions between the constituents. In this paper we investigate 
the relationship between four (previously disconnected) lattice results: i. spectral densities from 
MEM analysis of correlators; ii. static quark free energies F{R); iii. quasiparticle masses; iv. bulk 
thermodynamics p{T). We show a high degree of consistency among them not known before. The 
potentials V{R) derived from F{R) lead to large number of binary bound states, mostly colored, in 
gq, qq, gg, on top of the usual qq mesons. Using the Klein-Gordon equation and (ii-iii) we evaluate 
their binding energies and locate the zero binding endpoints on the phase diagram, which happen 
to agree with (i). We then estimate the contribution of all states to the bulk thermodynamics in 
agreement with (iv). We also address a number of theoreticall issues related with to the role of the 
quark/gluon spin in binding at large a,, although we do not yet include those in our estimates. Also 
the issue of the transport properties (viscosity, color conductivity) in this novel description of the 
QGP will be addressed elsewhere. 



I. INTRODUCTION 

A. New view of the Quark-gluon Plasma at 
not-too-high T 

From its inception two decades ago, the high-T phase 
of QCD commonly known as the Quark-Gluon Plasma 
(QGP) after [1], was described as a weakly interacting 
gas of "quasiparticles" (quarks and gluons). Indeed, at 
very high temperatures asymptotic freedom causes the 
electric coupling to be small and the QGP to be weakly 
interacting or perturbative *. At intermediate tempera- 
tures of few times the critical temperature Tc of imme- 
diate relevance to current experiments, there is new and 
growing evidence that the QGP is not weakly coupled. 
In a recent letter [3] (to be referred to below as SZl) 
we have proposed that in this region QCD seems to be 
close to a strongly coupled Coulomb regime, with an ef- 
fective coupling constant a « 0.5 — 1 and multiple bound 
states of quasiparticles. We have argued there (and will 
show it more quantitative below) that these bound states 
are very important for the thermodynamics of the QGP, 
responsible for roughly half of the pressure in RHIC do- 
main. We will not address transport properties of the 
QGP in this paper, although we expect the bound states 
to play an even more important role. 

To set the stage we first recall the chief arguments 
behind a scenario of strongly coupled QGP at such tem- 
peratures: i. a low viscosity argument; ii. recent lattice 



*The exception is long range color magnetism which remains 
nonperturbative [2]. 



findings of cc, qq bound states at T > T^, iii. the high 
pressure puzzle. 

Transport properties of QGP were so far stud- 
ied mostly perturbatively, in powers of the weak cou- 
pling. This approach predicted a large mean free path, 
Tlmfp ~ l/5^1n(l/g) :s> 1. Similar pQCD-inspired ideas 
have led to the pessimistic expectation that the Relativis- 
tic Heavy Ion Collider (RHIC) project in Brookhaven 
National Laboratory, would produce a firework of multi- 
ple mini-jets rather than QGP. However already the very 
first RHIC run, in the summer of 2000, has shown spec- 
tacular collective phenomena known as radial and elliptic 
flows. The spectra of about 99% of all kinds of secon- 
daries (except their high-pt tails) are accurately described 
by ideal hydrodynamics [4] . Further studies of partonic 
cascades [5] and viscosity corrections [6] have confirmed 
that one can only understand RHIC data by very low 
viscosity or large parton rescattering cross sections ex- 
ceeding pQCD predictions by large factors of about ~ 50 
or so. In short the QGP probed at RHIC is by far the 
perfect liquid known so far, with the smallest viscosity- 
to-entropy ratio ever, i.e. r]/s « .1 [6]. We note that 
from the theory standpoint ideal hydrodynamics, com- 
plemented by a "non-ideal" expansion in powers of the 
mean free path (the inverse powers of the rescattering 
cross section), is perhaps the oldest example of a strong 
coupling expansion. 

Naturally, these observations have increased our in- 
terest in other strongly interacting systems. Two such 
examples, discussed already in SZl are: i. trapped ul- 
tracold atoms driven to strong coupling via Feshbach 
resonances [7,8]; ii. J\f=4 supersymmetric gauge the- 
ory (CFT) recently studied via AdS/CFT correspon- 
dence [9,10]. In both cases, see the atomic experi- 
ments [11] for the first and the CFT viscosity calculation 
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in [12] for the second, strong coupling was found to lead 
to a hydrodynamic behavior, with a very small viscosity. 

The central point of the SZl paper was to provide at 
least a qualitative explanation to this small viscosity by 
relating it to multiple loosely bound binary states of quasi- 
particles, which should result in larger scattering lengths 
induced by low- lying resonances. We argued that at the 
zero binding points (indicated by lines in Fig. la and 
black dots in Fig. lb) those effects should be enhanced 
by a Breit-Wigner resonance, producing the cross sec- 
tion (modulo the obvious spin factors depending on the 
channel) 



a{k) 
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fc2 {E - Er)^ +T^JA 



(1) 



At E — Er ^ the in- and total widths are about equal 
and approximately cancel. The ensuing "unitarity lim- 
ited" scattering cross sections are large. This conjecture 
is nicely supported by the atomic experiments mentioned 
in [11], in which precisely the same mechanism was shown 
to ensure a hydrodynamical "elliptic flow" . 
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FIG. 1. Schematic position of several zero binding lines on 
the QCD phase diagram, from SZl. 

Bound states in the QGP phase. The earliest 
QGP signal (suggested by one of us [1]) was the disap- 
pearance of familiar hadronic states, especially the vec- 
tor ones - p, uj, 0, J /i}) mesons - directly observable via 
dilepton experiments. It was suggested later that even 
small-size and deeply-bound states of charmonium such 
as ?7c, J/i^ were expected to melt at T « Tc [13,14], so 
their absence was proposed to be a "golden signature" 
of the QGP. However, recent lattice works [15] using the 
Maximal Entropy Method (MEM) have found that char- 
monium states actually persists to at least T « 2 Tc, and 
there are similar evidences about mesonic bound states 
made of light quarks as well [16]. As we will show below, 
this is in good agreement with independent lattice stud- 
ies on the effective interaction potential between static 



sources in QGP . In SZl and now we argue that on top 
of familiar colorless hadronic states there should also be 
litter ally hundreds of colored binary bound states, and 
for the singlet gg pair (the most attractive channel) the 
region of binding should persist up to quite high temper- 
ature, i.e. r « 4Tr. 




1.0 1.5 2.0 2.5 3.0 3.5 4.0 
FIG. 2. The pressure normalized to that of a gas of mass- 
less and noninteracting quasiparticles, or Stephan-Boltzmann 
value, versus the temperature T /Tc, from the lattice calcu- 
lations by the Bielefeld group [18]. Different curves are for 
different numbers (masses) of the dynamical quarks. 

The "high pressure puzzle" stems from Fig. 2 
which shows a sample of lattice results for the QCD bulk 
pressure p(T), normalized to Stephan-Boltzmann nonin- 
teracting value for each theory. At T « 2 Tc the pressure 
nearly saturates to 0.8 of that for massless noninteracting 
gas, psB- 

This behaviour appears to be in contradiction with 
other lattice QCD works which found rather heavy quasi- 
particles, with masses (energies) e.g. of Mq^g/T ~ 3 at 
T Ri 1.5 — 2 Tf, [17]. How could such heavy quasiparticles 
account for large pressures [18]? Substantial elliptic flow 
effects at RHIC [4] points also to a large pressure in the 
prompt phase at RHIC or at T = (1 — 2) Tc, so we do 
know it is there. 

A similar discrepancy but now analytic and paramet- 
ric, was found for CFT at parametrically large cou- 
pling. In our second paper [19] we argued that in this 
case the matter cannot be made of quasiparticles, which 
are again too heavy parametrically, but rather by their 
(much lighter) binary composites. The QGP picture 
we advocate is thus just a beginning of this paramet- 
ric trend, with the intermediate running gauge coupling 
a, « (1/2-1). 



^Discrepancies with earlier results are mainly due to a con- 
fusion between a free and potential energy. 
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II. COLOR FORCES IN VARIOUS BINARY 
CHANNELS 

At T > Tc there is no color confinement, and so the 
interaction is a Coulomb-like at small distances, with a 
Debye-type screening [20] at large distances. Both lattice 
data as we will use below and results from the AdS/CFT 
correspondence agree that these features of the potential 
carry to the strong coupling regime. 

In this section we compare multiple binary colored 
channels ^, by using Casimir scaling for their relative 
strength. Its precise formulation can be made as follows. 
Let A,B he the color representations of either the quark 
or gluon constituents. If they are in a colored bound state 
with overall color representation D, then their color in- 
teraction is proportional to ^ 

cd={Xa-Xb)^ = 2{C{D)-C{A)-C{B)) , (2) 

where the C's are the pertinent Casimirs. For SU{3)c 
the Casimirs can be given in terms of the Dynkin index 
(mn) of the representation D, 

C{D) = m + n + ^{m^ + + mn) (3) 

In this section, we detail the color representations and 
the strength of the Coulomb interaction (2 in the bound 
states gg, qg and qq. 

gg : Two ghions yield the sum of irreducible color 
representations D = 8ig)8 = l©8s®8^®10s© 10 a ® 
27. In terms of the Dynkin indices, the same irreducible 
decomposition yields D = (00) © (11) © (11) © (03) © 
(30) © (22). Thus (2) reduces to 

CD= (a8-A8)^ = 2(C(£))-2C(8)) , (4) 

with Ci = —12, Cg = —6 (attractive) and Cio = (inac- 
tive) and C27 = 4 (repulsive). 

qg : A quark and a gluon yield the sum of irreducible 
color representations D = 3(g)8 = 3©65©15. In terms of 

the Dynkin indices, the same irreducible decomposition 
yields £> = (10) © (02) © (21). Thus (2) reduces to 

CB= (A3-A8)^ = 2(C(i))-4/3-3) , (5) 



*Multi-body bound states axe also allowed, but will not be 
discussed here. 

^This formula is analogous to the familiar SU(2) calculation 
of a relative spin projections in a state with some total spin 
J. 



with C3 = —6, Cg = —2 (attractive) and C15 = 2 (repul- 
sive). A similar decomposition applies to the conjugate 
representation qg with D = (01) © (20) © (12) 

qq : A quark and an antiquark yield the sum of ir- 
reducible color representations D = 3®3 = 1©8. In 

terms of the Dynkin indices, the same irreducible decom- 
position yields D = (00) © (11). Thus (2) reduces to 

Ci3 = (A3-A3)^ = 2(C(i?)-8/3) , (6) 

with ci = —16/3 (attractive) and cg = +2/3 (repulsive). 

qq . A two-quark state yield the sum of irreducible 
color representations _D = 3 3 = 3 © 6. In terms of 
the Dynkin indices, the same irreducible decomposition 
yields D = (01) © (20). Thus (2) reduces to 

CD = (As • A3) ^ = 2 iC{D) - 8/3) , (7) 

with Cg = —8/3 (attractive) and Ce = -1-4/3 (repulsive). 

Using a singlet qq as a standard benchmark (the only 
one studied extensively on the lattice), one can sum- 
marise the list of all attractive channels in Table I, in- 
dicating the relative strength of the Coulomb potential 
in a given color channel and the number of states. Even 
without the excited states to be discussed below, there is 
a total of 481 channels for two flavors (41 colorless), and 
749 states (81 colorless) for three flavors. 

III. BOUND STATES IN STRONG COULOMB 
FIELD 

Before we analyze the relativistic two-body bound 
states for quarks and gluons, we go over the results for the 
simpler problems involving cither a spin 0, spin 1/2 and 
spin 1 particle moving in a strong (color) Coulomb field, 
where the effect of color is treated as a Casimir rescal- 
ing of the Coulomb charge. Precession in color space can 
be treated but will be ignored through an "Abelianiza- 
tion" of the external field. The results for spin and 1 /2 
are known from 1928 [21]. They are presented for com- 
pleteness since they streamline our analysis for spin 1. A 
canonical application of the latter is that of a W boson 
bound to heavy Coulomb center. 



channel 


rep. 


charge factor 


no. of states 


99 


1 


9/4 


9s 


99 


8 


9/8 


9, * 16 


99 + 99 


3 


9/8 


3c * 6s * 2 * Af/ 


19 + 99 


6 


3/8 


6c * 6s * 2 * Nf 


m 


1 


1 




99 + Q9 


3 


1/2 


4s * 3c * 2 * ATj 



TABLE I. Binary attractive channels discussed in this 
work, the subscripts s,c,f mean spin, color and flavor, Nf = 3 
is the number of relevant flavors. 
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A. Spin 



B. Spin i 



For a scalar particle the Klein-Gordon equation 



{E - gVy -m^+d i> = 



(8) 



should be used. This equation was analyzed for a 
Coulomb potential** gV = —a/r in [21] with the energy 
spectrum ^ 



E{nr, I) 



1/2 



m 



1 + 



K + + 1/2)2 -a2)2 



= 1 (9) 



Taking the lowest level to be = 1, Z = as an example, 

one finds that a = 1/2 is a critical value for this equation. 
Although the binding is at this point finite and even not 
large, E{l,0)/m = a/4/5, something new is obviously 
happening at this critical coupling because the square 
root (in the denominator) goes complex. 

What happens is that the particle starts falling towards 
the center. Indeed, ignoring at small r all terms except 
the V"^ term one finds that the radial equation is 
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(10) 



which at small r has a general solution 



R = Ar"* + Br"- , s± = -1/2 ± ^/T/4~~a2 (11) 

that for a 1/2 is just l/r^/^. At the critical coupling 
both solutions have the same (singular) behavior at small 
r. For a > 1/2 the falling starts, as one sees from the 
complex (oscillating) solutions. 

In the CFT theory with fixed (non-running) coupling 
constant, nothing can prevent the particle from failing 
to arbitrary small r as soon asa>l/2M_ in con- 
trast, in QCD the coupling runs, a{r) w l/ln(l/rAQCD)5 
so that at small enough distances the coupling gets less 
than critical and the falling stops. In [22] a crude model 
of a regularized Coulomb field was used, producing the 
same effect. Our arguments show that asymptotic free- 
dom would be in principle enough. However, the wave 
function at the origin is changing dramatically at a w 1/2 
and in view of that we performed an additional study of 
the Klein-Gordon problem with a Coulomb-|-quasi-local 
potential in the Appendix. 

The falling onto the center happens for any spin of 
the particle, only the value for the critical coupling is 
different. We now proceed to show that. 



**Wc use in this section a QED-like notations in which a = 
. The standard QCD notations for fundamental charges 
should read a = {A/3)g^/4n = {'i/3)a,. 

^hn [19] a WKB analysis was used with apologies to [21]. 

•'"'"In fact, this is why the dual string description has a black 
hole. One of us (ES) thanks Daniel Kabat for pointing this 
to him. 



The squared Dirac equation for a massive spinor * 
reads 



(12) 



with a^^ = i[7ju,7y]/2 and F'^" the external background 
field. In the chiral basis, (12) simplifies 

((□+m2) + 25S-(BTiE))*±=0, (13) 

where the spin operator is Lie-algebra valued, 

[S°,S''] =ir^^S^ (14) 

The magnetic contribution in (12) is standard. The elec- 
tric contribution is complex and is reminiscent of a 
Bohm-Aharanov effect. 

The rclativistic Coulomb problem stemming from (12) 
for each of the two stationary spin components reads 




the solutions to (15) are naturally sought in terms of 
spinor spherical harmonics [23] 



J L k \ etY^-'' , (16) 



-M {M -b) b I ^ 



with the bracket a conventional Clebsch-Gordon coeffi- 
cient restricting the values of L and a, & = ±1/2. In this 
representation the spinors arc just = (J"**, and (16) is 
an eigenstate of J^, J3 and with eigenvalues J{J +1), 
~M and L{L + 1) respectively. In the basis (16) the spin 
operator S • f is off-diagonal 



(17) 



Using the spinor spherical harmonics (16) in (15) 
through the expansion 



(18) 



for fixed JM and (17) we obtain the 2x2 matrix equation 
for the radial function R 

f^^9\2 2^cP 2d J{J +1) + C \ „ 

[E H y -m^ + + ^ 5 I R = 

r ar^ r ar 



Recall that the squared Dirac operator is hermitean. 



with 



-W 4 

The eigenvalues of C are 



(20) 



(21) 



In terms of (21) the cigcnvahie equation (19) becomes 
diagonal for the rotated R. Defining e = \/m?- — E'^ and 
r = a;/2e and R = 2eU/\/S yield (19) in the diagonal 
basis in the form 



with 



a;U" + U'+(a-^-^) U 
4 4a;/ 



a = g — 



(22) 



(23) 



The equation (22) is a standard hypergeometric equation. 
The bound state solutions with \E\ < m are Laguerre 
polynomials for 



a 



1 /? 
-+2+2 



(24) 



which is the quantization condition with integer n refer- 
ring to the nodal number of the wavefunction as opposed 
to the radial quantum number n^. used above. Unwinding 
(24) in terms of (23) yields the spectra for the relativistic 
spin 1/2 particles/antiparticles 




= ±1 



(25) 



with (23) given explicitly by 
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(26) 



after using (21). This result is in agreement with the 
original result established in [21]. Note that for scalars 
it agrees with the semiclassical results we used earlier. 



C. Spin 1 



□ 5^. - (1 - |) V. - 2igF^^ ) a"^ = (27) 



which simplifies for Feynman gauge ^ = 1 to 



(28) 



with y ij,af^ = 0. The gluon in (28) has two physical 
polarizations, the longitudinal and time-like ones being 

gauge artifacts. Generically, we can decompose the gluon 
along its polarizations a'* = e(^*", and rewrite (28) in 
the form 



(29) 



with «S^j^/4 



T 



This result is reminiscent of the 



massless spin i equation (12) if we were to interpret 
as the spin-operator of the gluon in the polarization space 

For massive gluons the pertinent equation is a variant 

of the Proca equation extensively used in the litterature. 
Here instead, we proceed by analogy with the spin 1/2 
case. In the polarization space the equation of motion 
for spin 1 is just (the g- factor is now 1 instead of 2) 



((□ + m^) + 5 S • (B - iE)) * = 



(30) 



which is readily checked from the path integral approach 

in the first quantization (see also below). Note that 
(30) is a 3 X 3 matrix equation for one longitudinal 
and two transverse polarizations. For spin 1 gluons, 
gAafa _ jgAab an external electric gluon field, (30) 
reduces to 



* + 5 (E X *) = 



(31) 



which shows that the electric field causes the polarization 
to precess in the relativistic equation. To solve (30) in 
a Coulomb field E = —gfjr^ we use exactly the same 
method discussed for spin i except for the use of vector 
instead of spinor spherical harmonics, 



* JML 



J 

-M 



L 

{M-h) 



M-b 



(32) 



with a,b = 0, ±1. In this representation the three polar- 
izations are chosen real with again = 5"^. A rerun of 
the precedent arguments show that (19) holds for spin 1 
in a 3 X 3 matrix form with the substitution 



/l 



C = 



-ig 



(2J+1) 
/ J+i 

V 2J+1 





-ig- 



-tg 



j+i 

2J+1 







2J+1 

i + (2j+i)y 



(33) 



The same analysis performed above can be carried out 
for the gauge-independent part of the wavefunction. In- 
deed, let us first consider a massless gluon in an arbitrary 
covariant background field gauge. The equation of mo- 
tion is standard and reads 



***Thc same equation (29) follows from a path-integral de- 
scription of a quantum mechanical evolution of a massless 
spin-1 particle in which is the covariantized spin. 
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for J ^ 0. The case J = is special since C = 
diag(0, 0, 2). The eigenvalues of C are solution to a cubic 
(Cardano) equation 

+ + A(l - (2 J + 1)2 + + 2/ = . (34) 

The solutions are all real since the polynomial discrimant 
of (34) is negative [24]. This is expected since the gluon 
evolution operator is hermitean. The explicit solutions 
to (34) are [24] 



-Qcos 



9 + 2-Ka 



(35) 



with cos^ = R/ \J —Q^ and 

Q = i(l-(2J+l)2+/)-l 

i?=i(l-(2J+l)2+5^)-/-A 



(36) 



The corresponding spectrum for spin 1 particle is again 
of the type 



^ 1+ 

m\ (n+i + |)_ 



1/2 



±1 



with 



(37) 



(38) 



The case J = is special and yields A = (0, 0, 2). (Thus 
one can see that the minimal critical coupling is q = 1/2, 
as for spinless case.) Much like the quarks and scalars, 
the gluons fall onto the center at a critical coupling which 
is now set by the branch point of not only (38) but also 
(35). 



IV. A TWO-BODY BOUND STATE PROBLEM 
A. Generalities 

Non-rclativistically, the problem of two bodies (e.g. 
positronium) is readily reduced to a single one body prob- 
lem (hydrogen atom) by a simple adjustment of the par- 
ticle mass. Relativistically this is more subtle. 

Let us start with two spinless particles, obeying two 
separate Klein Gordon equations. Making use of the fact 
that in the CM the two momenta are equal and opposite 
p — p{l) = — 1?(2) one can eliminate the relative energy 
Po(l) -Po(2) for the total energy E = po{l) +Po(2). The 
resulting relations 

K^-f = 

«2 = {E/2r - (m? + mi)/2 + ^^^^i^^ (39) 



can be used as the KG equation for two particles with dif- 
ferent masses. The KG simplifies for equal masses, when 
the last term drops out. Fortunately this is is approx- 
imately true for our problem, since for T in the region 
of interest the quasiparticle mass difference is relatively 
small even for qg states. 

The quantization follows from the canonical substitu- 
tion E — > —idt — V, p ^ —id — gA. The magnetic 
effects through A, are comparable to the electric effects 
for w pa 1. 




FIG. 3. ClEissical motion of two mutually attracting 
charges. If the field propagation is instantaneous, a partner 

of a particle P is at the "antipodc" point A. However if parti- 
cles move relativistically and the electromagnetic field travels 
at light speeds, the particle P sees a field from the earlier 
position B (drawn for counterclockwise rotation). 

The next problem is to understand what exactly is the 
4 potential {V, A) in this equation. That is of course the 
field at one particle (say P in Fig. 3) due to the other 
one. Non-rclativistically the particle speed is negligible 
compared to the speed of light, so one can safely place 
the other particle at the opposite point (say A in Fig. 3). 

Including retardation, one classically expects the field 
to emanate from B instead of A. The retarded position 
B is simply determined by the condition that the travel 
time from B to A equals the time it takes light from B to 
P, that is {AB)/v — {PB)/c. In the ultrarelativistic case 
V = c and one gets a simple equation for the max;imal 
retardation angle a, 



a = \/2 + 2 



cos a 



(40) 



with a root at a = 1.48. The retardation angle is about 
85° which is rather large. 

The classical description is oversimplified, and in fact 
quantum theory allows for field propagation with any 
speed, not just c. It demands a convolution of the 
path (current) with the quantum propagator of a pho- 
ton (gluon). At this point, it may appear that all hope 
to keep a potential model is lost and a retreat to a full 
quantum field theory treatment is inevitable. 

This is indeed the case for intermediate coupling a « 1, 
but when it gets stronger one can again claim some virtue 
in a potential-based approach. Our argument presented 
in [19] resulted in the conclusion that the effective photon 
(gluon) speed is not the usual speed of light but larger, 
by a parameter ^\/A » 1 (where as usual A = g'^Nc). 



6 



Thus the field is stiU dominated by the emission at the 
"antipode" that is point A instead of B. 

The magnetic effects are the usual current-current in- 
teractions, present for spinless particles as well, plus 
those induced by (gluo) magnetic moments related with 
particle spins, plus their combination (spin-orbit). The 
spin effects were argued to be small, see [22,19]. 

For an extensive review of the known results on how 
one can reduce a 2-body relativistic Dirac problem to 
that of a potential problem we refer to [25]. Here we just 
note as in [19] that in the QGP the quark mass is a "chiral 
mass" , so the derivation of the effective single-body Dirac 
equation in this case would be a priori different from the 
one discussed in [25]. 



B. Relativistic Two-Body Bound States 



(46) 



which is the same as 



2 



(47) 



This form is standard, with the Coulomb (Gauss law) 
and the current-current (Ampere's law) interaction. For 
two particles the spin effects are encoded in the spinors 
^ along with the particle-antiparticlc content. They 
may be unravelled non-relativistically using a Foldy- 
Woutuhysen transformation. 



In this section we illustrate the derivation of the rela- 
tivistic two-body bound states between qq, gg and qg in 
QCD by considering the simpler problem in relativistic 
QED in both the first and second quantized form. The 
generalization to QCD is staightforward in the canonical 
quantization approach i.e. = gauge. 



1. Second Quantization 

Consider two massive relativistic electrons with Dirac 
spinors \l/a and a = 1, 2 



£ = 5^ (i)(7 - m„) - -F^^F^^ . 



(41) 



Canonical quantization of (41) yields the hamiltonian 
density 

n = nF + ^(E2 + B2) - J. 1+ (Jo - a • E) Ao (42) 
with Hf the free fermion hamiltonian density and 

J" = 5 ^ *a7'^*a . (43) 
a=l,2 

The coefficient of Aq (constraint) is Gauss law. The latter 
is resolved in terms of the Coulomb potential 



(44) 



Ampere's law follows from the equation of motion for the 
transverse part of the electric field 



E7 



H , E J" 



d'^Ar-J. 



(45) 



2. First Quantization 

Perhaps a more transparent way to address the spin 

effects in the presence of gauge fields when the particle- 
antiparticle content is not dominant, is to use the first 
quantized form of the same problem. For that it is best 
to choose the einbein formulation [26] in the rest frame. 
For two gauge coupled particles it reads 

c = -^{i-ii)-^^+aA°-aS,-A. 
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(48) 



where a = 1,2 is summed over throughout unless in- 
dicated otherwise, and Aa = A{xa) and similarly for 
the fields. The einbeins are denoted by Va- They will 
be considered as Lagrange multipliers at the end and 
fixed by minimizing the energy. The =F refers to particle- 
antiparticlc. 

We note that while the potentials couple canonically 
to the currents, the spin couples canonically to the mag- 
netic field but has an imaginary coupling to the electric 
field in Minkowski space, a situation reminiscent of Berry 
phases. Indeed, canonical quantization of (48) yields the 
momenta (only the upper sign is retained from now on) 



Pa = VaXa- gAa 



(49) 



The canonical momentum 11 has a complex shift due to 
the spin of the particle in contrast to the second quan- 
tized analysis. By insisting that (49) are canonical, we 



For stationary (bound) states At = (1/3^) J and the 
hamiltonian density (42) simplifies 



tt^The particle content of the Lagrangian lives only on the 
time axis. 
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conclude that the energy spectrum is shifted by the Berry 
phase. 

The canonical Hamiltonian following from (48) after 
resolving Gauss law reads 



1 



Va 



s ^ ^ 

-g-7J- ■ (Ba - il^La - fflTa) 



S„ -8,5(0) (50) 



with 



' Q2 



(51) 



with the summation over a shown explicitly. 

For stationary states the equation of motions can be 
used to solve for At and therefore for IIt and B through 
Ampere's and Lenz's law. In particular, 



At{x) = - 



1 



Vad^ 
1 



{Pa+ gATa)&{x ~ Xa) 

g{^a ^ d) 5{x - xa) . (52) 



Noting that the momenta scales with the einbeins as 
Hi « V° and that w B w it follows that 

for stationary states the hamiltonian simplifies to order 
0{l/V^), i.e. 



7i ; 



1 



2igSa ■ Ilia 



1 



(53) 



The expansion in is justified in the non-relativistic 

limit since V ~ m and the ultrarelativistic limit since 
V « (here 7 is the Lorentz contraction factor). In- 
deed, the extrema in V of the total hamiltonian H asso- 
ciated to (53) are complex and read 



Va = \/pI 



ml + 2igSa ■ ULa 



(54) 



in particular V ~ 7m as asserted. For a pair of identical 
particles mi = m.2 in their center of mass frame pi — 
—P2 = P, (53) and (54) yield 



4 1 X12, 



- 2i.g2(Si - S2) • d— (55) 
X12 



after absorbing the Coulomb self-energy corrections in 
the masse m . (55) gives a spectrum analogous to the 
one described in the background field section except for 
the fact that now each of the two particles carry its own 
spin in the center of mass frame. Spin-orbit and spin- 
spin effects are obtained by carrying the expansion a step 
further in 1/V in H. 



V. BOUND STATES IN STATIC EFFECTIVE 
POTENTIALS 

Studies of effective potentials in lattice QCD have a 
long history. Their T — version was first obtained in 
the classic paper by M. Creutz who first found confine- 
ment on the lattice in 1979. The first finite-T results 
have shown Debye screening, in agreement with theo- 
retical expectations [1]. These static potentials lead to 
early conclusions [13,14] that all states, even the lowest 
cc states and r]c, J/ip, melt at T « Tc- As we already 
mentioned in the introduction, these conclusions contra- 
dict the recent MEM analysis of the correlators which 
indicate that charmonium states stubbornly persist till 
about 2Tc. 

On theoretical grounds, it has been repeatedly argued 
(see e.g. [27]) that close to Tc the Debye mass is low 
enough to allow the color charge to run to rather large 
values. If so, the binding of many states must occur, as 
it was shown in our letter [3]. 

Recently in a number of publications, the Bielefeld 
group had reanalyzed the effective static potentials. One 
of the key point missed in previous works is that the ex- 
pectations values obtained from thermal average are by 
definition the free energy F, related with the energy E by 
the standard thermodynamical relation 



E = F -TS 



(56) 



So in order to understand the potential energy one has to 
remove the entropy part first. The subtraction results in 
in much deeper potentials, which (as we will show below) 
readily bind heavy (and light) quarks. 




1 2 3 4 5 6 

FIG. 4. Heavy quark free energies in the singlet channel for 
2-flavors of dynamical quarks at a quark mass of m/T = 0.40 
on 16"' X 4 lattices renormalized to the zero-T potential ob- 
tained from [28]. 

A set of potentials obtained by the Bielefeld group is 
shown in Fig. 4. The strength of the effective interaction 
can be characterize by a combination, called a screening 
function, 
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^i(r,T) = --r(Fi(r,T) 



i^i(oo,r)) 



(57) 



where the subscript 1 refers to the color singlet channel 
and the 3/4 removes the Casimir for qq representation, 
so that Si is in a way just an effective gauge coupling 
as . A sample of these effective gauge couplings are plot- 
ted in Fig. 5. The logarithmic plot shows exponential 
decrease at large r, while the non-logarithmic plot shows 
a decrease toward small r due to asymptotic freedom. 
The maximum at rT « 1/2 indicates that the effective 
Coulomb coupling at is a^s ~ 4/35*1 ~ 1/2, right in 
the ballpark used in [22]. One should also note that the 
strength is even larger below Tc, where it is related with 
confinement at the string breaking point. Finally, we 
note the surprisingly monotonous behavior through the 
phase transition point, due to the fact that at T > Tc 
the static charge continues to be screened by a single 
light quark, like in a heavy-light (B-like) meson below 



1.4 



1.2 

1 

0.8 
0.6 
0.4 
0.2 




Si(r,T) 



X ^ tr 



TAT.=6.75 
T/r=0.82 
T/T.=0.91 
T/T^=0.97 
T/T.=1.01 

Tn"c=1 .07 
T/Tt.=1.12 
T/T.=1.16 
T/Tp=1.24 



Hit, 



0.2 0.4 0.6 0.8 1 1.2 1.4 

FIG. 5. The color singlet screening function, from [28]. 

The effective quark mass, or a constant value of the 
potential, was subtracted out: it plays an important role 
in what follows. 

We have parameterized these Bielefeld data by the fol- 
lowing expression (here and below all dimensions are set 
up by Tc, e.g. T means T/Tc and r is rTc) 



Fi{T,r) = 1.5 -1.1{T-1) 



1/2 



3r ln(l/r + 3T) 



(58) 



and then extracted from it the potential energy using 
Ei=Fi- TdFi/dT. 

Furthermore, an appropriate normalization of the po- 
tential to be used for the bound state analysis is its decay 
at infinity We use V{T,r) = Ei{T,r) - Ei{T,r ^ oo). 
(Non-relativistically, the constant part can be added to 
a mass.) The resulting potentials are plotted in Fig. 6 

Our next next step after extracting the lattice poten- 
tials 6 is to use them in a Schrodinger (or appropriate 
relativistic) equation and solve for the bound states. A 



sample of the results obtained using the Klein-Gordon 
equation with this potential is plotted in Fig. 7. We used 
a charm quark mass of 1.5 GeV and an effective gluon 
mass of 0.6 GeV: the results are not very sensitive to it. 




n — I I I I I — I I I I — I I I I I 
0.2 0.4 0.6 



I I I I 
. 



FIG. 6. The static potentials V(T,r) as a function of 
the distance rTc. The values of temperature used are 
T = 1., 1.2, 1.4, 2., 4., 6., lO.Tc, from right to left. 

One can see that charmonium remains bound to about 
T = 2.7Tc. It is in fact completely consistent with lattice 
observations [15] using MEM. The fact that the state is 
traced only to T < 2Tc is completely understandable, as 
for 2 < T/Tc < 2.7 it is so weakly bound that the size 
of the state may exceed the size of the lattice and could 
not possible be seen. Note that at all T the charmonium 
binding remains rather small, and so the nonrelativistic 
treatment of charmonium would be completely justified. 
A similar conclusion would be reached for light qq pairs. 

This is not the case for the gluonic singlet gg state, 
which has a color charge larger by the ratio 9/4. We 
found that the same potential in this case leads to much 
larger binding, reaching up to 40 percent of the total 
mass at T = Tc. There is no question that the relativistic 
treatment is indeed needed here. 

We have also looked for Z = 1 states in this potential, 
which we only found for the most attractive singlet gg 
state. Those reache their zero binding at T = 1.0805 Tg. 
Another next-shell state is an Z = 0, n = 2 state, which 
exists till T < 1.205 Tc. (Those states are not included 
in the calculation of the thcrmodynamical quantities at 
the end of the paper.) 
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VI. RELATIVISTIC CORRECTIONS AND 
INSTANTON MOLECULES 

A. Relativistic corrections 

The lattice potentials used above were evaluated for 
static charges without spin, therefore it docs not include 
effects proportional to particle velocities and/or spins. 
Although we used a relativistic KG equation, we see that 
qq (and even the most attractive color singlet gg state) 
arc not bound nearly enough as to become massless. On 
the other hand, we know that the lowest qq states, the 
pion-sigma multiplet, must do so at T = Tc. It means 
that something is missing in the interaction. We will 
discuss those missing effects in the next section. 

The first relativistic effect, already discussed in 
Ref. [22] , is the velocity- velocity force due to magnetic in- 
teraction (Ampere's law in the classical treatment). This 
corresponds to a substitution of the effective Coulomb 
coupling by 

a ^ a(l - P1P2) (59) 

where (3i are velocities (in units of c) of both charges. In 
the center of mass the velocities are always opposite, so 
the the effective coupling always increases. 

We have estimated the mean velocity squared by using 
the equation itself 

<«'>=/*^»(^-(e^)1 (60) 

where the wave function is appropriately normalized. 
This yields < >= .12 for light quarks at T/Tc = 1.05. 
If we plug this correction back into the potential, as- 
suming it scales as (1+ < >) as a whole we get 
larger binding (squares in Fig. 7). With this relativistic 
correction included, at T w Tc the light quarks become 
relativistic and about as bound as the charmed quarks, 
with mean velocity of about 1/3. 

Although we do not yet see how very relativistic mo- 
tion may come about, we will do so in the next section. 
In anticipation of that, let us show here what can be 
a maximal effect of the the relativistic correction under 
consideration. If the particle velocity becomes the speed 
of light, a correction imdcr consideration effectively dou- 
bles the coupling, putting it at T r:! to be aeff ~ 1. If 
we simply double the effective potential as a whole, the 



we argued above, Ampere's current-current interac- 
tion is times the Coulomb charge-charge interaction, but 
in order for the whole screened static potential to scale the 
same way, other parameters (such as the electric and magnetic 
screening lengths) should coincide, which strictly speaking is 
not the case. 



binding increases significantly. The result is shown by 
a single square in Fig. 7(a), and ai T = Tc the binding 
reaches about 1/4 of the total mass. Even larger effect 
is seen in the particle density at the origin. As shown in 
Fig. 7(b) it increases by about an order of magnitude. 
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FIG. 7. (a) The energy of the bound state E (in units of 
the total mass 2M) versus the temperature T/Tc using the lat- 
tice effective potential "(/^(T, r), for charmoiiium (crosses and 
dashed line), singlet light quarks qq (solid line) and gg (solid 
line with circles). A set of squares show the relativistic cor- 
rection to light quark, a single square at T = 1.05Tc is for qq 
with twice the coupling, which is the maximal possible rela- 
tivistic correction, (b) The density at the origin |i/)(0)p/r^ 
of the bound states versus the temperature T/Tc in the lat- 
tice effective potential V(T, r), for singlet light quarks qq and 
gg (upper and lower lines with circles). A single square at 
T = l.OSTc is for qq with twice the coupling. 
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B. Interaction induced by the instanton molecules 

We have seen in the preceding section that relativistic 
effects proportional to velocities make all states signifi- 
cantly more bound and dense at T ^ Tci but still these 
effects are too weak to bring the total energy of the low- 
est qq to zero, as is required for sigma mesons to trigger 
a phase transition at T = Tc coming from above on the 
temperature axis. 

As discussed in BLRS paper [22], the missing element 
is a quasi-local interaction induced by the instanton- 
antiinstanton molecules, the lowest clusters of zero total 
topological charge allowed in the chirally restored phase 
(see [29] for a review). 

In this paper we will not try to estimate the cou- 
pling from first principle but instead adopt a purely phe- 
nomenological approach. For the local 4-fermion inter- 
action with the coupling constant G the energy shift is 
given simply by 

SE = -G\i>{0)f . (61) 

Here we will tune the magnitude of the effective 4-fcrmion 
coupling G^q so that the pion-sigma multiplet gets mass- 
less (and then tachyonic) exactly at T = Tc- 

With \SE\ = 0.92Mq « 1.3GeV and |V(0)|2 ~ 350T^, 
after the relativistic correction is included, one needs 
G « l.SGeV^ which is in the expected ballpark If 
the relativistic correction would not be there, the value 
of |V'(0)P would be an order of magnitude smaller, see 
Fig. 7(b). 

The remaining problem is what happens in the glueball 
{gg singlet) channel, where we expect the interaction with 
instantons to be even stronger than with quarks. Indeed, 
the instantons arc classical objects made of gluon field, 
thus their interaction with gluons would be proportional 
to the action 0{S w Sn^/g"^) which is expected to be 
about 10 times stronger than 't Hooft's interaction with 
quarks. If this is the case, the s-wave glueball state gets 
tachyonic already at some finite T above Tc- 

VII. CONTRIBUTION OF THE BOUND STATES 
TO THERMODYNAMICAL QUANTITIES 

A. The contribution of the quasiparticles 

Let us start with the papers by Pcshier et al [30] and 
Levai and Heinz [31] who used a simple quasiparticle gas 
model, deducing what properties the q,g quasiparticles 
should have in order to reproduce the lattice data of the 



pressure * p{T). Assuming the usual dispersion relation 
o;^ = p2 -I- M{T)-^, one has to deal with the T-dcpendent 
of the masses. For example, for pure glue it is possible 
to reproduce p{T) by assuming Mg(T). The qualitative 
behavior found in [31] at high T is about linear and rising, 
as expected perturbatively [20]. For T = (1.3 - 3)Tc it 
is about constant, with a rise towards Tc, an indication 
of the onset of confinement. 

Although these features arc qualitatively consistent 
with the direct measurements [17], they are not quan- 
titatively. For Nf = 2 the expected masses at T = 1.5 Tc 
(close to their minimum) from [31] are 

Mg w 420 MeV Mq w 300 MeV . (62) 

Howcvcir direct studies by Karsch and Laerman found 
heavier ones 

Mg « 540 MeV « 620 MeV . (63) 

If the reader is not impressed by this difference, let us 
mention that the corresponding Boltzmann factors for 
quarks are exp{-Mg/T) = 0.28 for LH and only 0.075 
for KL values. This means that the QGP quasiparticles 
at such T are too heavy to reproduce the global thermo- 
dynamical obscrvablcs. 

If this numerical example is not convincing, let us go 
to A/'=4 supersymmetric Yang-Mills theory at finite tem- 
perature, for which a parametric statement can be made. 
At strong coupling A = g^Nc » 1 the quasiparticle 
masses are [32] m « VXT T, and thus the corre- 
sponding Boltzmann suppression is about exp(— #\/A). 
In our work [19] we suggested that in this limit the mat- 
ter is made entirely of binary bound states with masses 
m ~ T. We have shown that such light and highly rela- 
tivistic bound states exist at any coupling, balanced by 
high angular momentum I « y/X. Furthermore, we have 
found that the density of such states remains constant 
at arbitrarily large coupling, although the energy of each 
individual state and even its existence depend on A. So, 
in this theory a transition from weak to strong coupling 
basically implies a smooth transition from a gas of quasi- 
particles to a gas of "dimers" . In QCD going from high 
T towards Tc the coupling changes from weak to medium 
strong values = 0.5-1 with A = 15-30, and so one nat- 
urally can see only a half of such phenomenon, with the 
contribution of the bound states becoming comparable 
to that of the quasiparticles. 

B. Parameterizations for masses of the bound states 

As we indicated above, the pressure problem can be 
solved when one accounts for the additional contribu- 



^^^The estimated value is a factor of 2 smaller than used 
in [22], which was derived by continuity from NJL-like fits 

from T < Te, and somewhat overestimated. ^^^^"^ t^*^'^"^*! observables follow from this function, 

we will not discuss them. 



11 



tion of the binary bound states. We will do so in a 
slightly schematic way, combining them into large blocks 
of states, rather than going over all attractive channels 
of Table 1 one at a time. 

The "pion multiplet" (plus other chiral and Ua{^) 
partners, e.g. (J,ri,S for Nj — 2) carry 2Nj states, which 
are to turn massless at Tc. Using (61) for the pion bind- 
ing and simple parameterization of the T-dependence of 
the wave function at the origin as shown in Fig. 7(b), we 
arrive at the following parametrization of their effective 
mass at temperature T 



R'U ~ 10 T, (1 - exp(-3(r - T,)/T,)) 



(64) 



which is set to vanish at T = Tc- 

For the "rho multiplet" (vectors and axial mesons) 
with 6Nj qq states we use the same expression, but with a 
suppression factor of 0.7 in front of the exponent (see [22] 
for details on why instanton molecules are somewhat less 
effective in vector channels). 

We have further ignored the (most bound) gg singlets 
and concentrated on much more numerous colored at- 
tractive channels. Ignoring the differences between them 
for simplicity, we lump them altogether with a mass pa- 
rameterized as 

Meoiored ~ 11-5 T, ((T/3T,)°-5 + Q.ITJ{T - T,)) (65) 

where the second term enforces their disappearance at 
the critical point. The corresponding curves for the 
masses of bound states are shown in Fig. 8a. 



C. The QGP pressure 

In analogy with the so called "resonance gas" at T < 
Tc, the contributions of all bound states for T > Tc can 
be simply added to the statistical sum, as independent 
particles. However, this is only true for sufficiently well 
bound states. The region near the zero binding point 
needs special attention. Indeed, in this region the bound 
state becomes virtual and leaves the physical part of the 
complex energy plane: one may naturally expects that its 
contribution to the partition function is also going away. 

Let us see how this happens, restricting ourselves to 
the s-wave (1=0) scattering*. As it is well known (see e.g. 
Landau-Lifshits QM, section 131) a general amplitude 
for a system with a shallow level can be written at small 
scattering momentum k as 



1 



—Kg + rok^ /2 — ik 



(66) 



*The case of Z > is different, as the centrifugal barrier 
allows for quasistationary states and resonances to exist. 



where kq is related to the level position and rg is the 
so called interaction range. The subscript refers to the 
I — partial wave. The cross section is 



47r 



{ko - rofc2/2)2 + F 



(67) 



and if the range term is further ignored one gets he fa- 
miliar a = At: /•m{E + |e|) form which does not care for 
the sign of the binding energy e = Kg/m. 
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FIG. 8. (a) The lines show^ twice the effective masses for 
quarks and gluons versus temperature T /T^. Note that for 
T < 3Tc Mq > Mg. Circles and squares indicate estimated 
masses of the pion-like and rho-like qq bound states, while the 
crosses stand for all colored states, (b) Pressure (in units of 
that for a gas of massless and noninteracting quasiparticles) 
versus the temperature T/Tc- The crosses correspond to the 
Nf — 2 lattice results, from Fig. 2. The lower solid curve is the 
contribution of unbound quasiparticles, the upper includes 
also that of all bound states. Squares are for the pion-like 
and rho-like qq bound states combined, and circles for all the 
colored bound states. 
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As kq goes through zero and changes sign, the scatter- 
ing length f{k = 0) = — I/kq jump from —oo to oo. As 
one can see from these expressions, a bound level close 
to zero generates significant repulsive interaction of the 
quasiparticles (with positive energies). As we will see 
shortly, an account for such repulsion reduces the contri- 
bution of the bound state to the partition function. 

Let us follow the well known Beth-Uhlenbeck expres- 
sion for non-ideal gases (see e.g. Landau-Lifshitz, Statis- 
tical mechanics, section 77) 



1 r° 



^e-^^'-^dk (68) 
dk 



where the first sum runs over all bound states with bind- 
ing energies e„ and the second over the scattering states. 
As the simplest example, let us consider the zero binding 
point Kq = 0, for which the expression for the scattering 
phase can be simplified to 



exp(2i^o(fc)) = 



rok/2 + i 
rofc/2 — i 



(69) 



and assuming that the temperature T is high enough to 
ignore the Boltzmann factor in the integral one gets 

Zint ~ 1 + -((5o(oo) - ,5o(0)) = 1/2 . (70) 

TT 

Thus, at the zero binding point the repulsion reduces the 
contribution of the bound states to a half its value. As 
the virtual level moves away from zero, the contribution 
decreases further'''. 

We will use a simple Fermi-like function to enforce this 
disappearance of the level from the statistical sum, multi- 
plying the level contribution by an additional "reduction 
factor" 



R{T) 



1 



l + exp[C(T-T,.6.)] 



(71) 



which reduces the contribution at the zero binding point 
Tz.b. by 1/2 and eliminates it at higher T. We will use a 
parameter C = 2/Tc. 

Taking it all together, one finally gets the pressure 
shown in Fig. 8b. One can see that the pions and rhos 
peak at T^, but become unimportant for higher T. The 



^Note that this is different from a contribution of a real res- 
onance, which generates a Breit-Wigner cross section with a 
maximum, and for a narrow resonance the same contribution 
as for the bound state. 

■'■Recall that the pions gets massless in the chiral Mmit. Prom 
the current description it looks like there is a disagreement 
with the lattice data at T = Tc, but one should recall that 
the latters are for medium heavy quarks. 



colored states are too heavy near Tc, but then their con- 
tribution to the pressure becomes comparable to that 
of the quasiparticle gas. There arc hundreds of them, 
bringing quite a large statistical weight, which is however 
tamed by a large mass and consequently small Boltzmann 
factors. As a result, the total contribution follows rea- 
sonably well the lattice data points. 



VIII. CONCLUSIONS AND DISCUSSION 

In this paper we have addressed a number of issues 
related with the bound states in the QGP phase at not 
too high temperature. We have catalogued all attractive 
binary channels, with proper color factors and multiplici- 
ties, see Table 1. We have also presented a unified frame- 
work for analysing two body relativistic bound states of 
arbitrary spin and mass using first quantization argu- 
ments. 

As already mentioned in the Abstract, our main task 
was to check the internal consistency of four (previously 
disconnected) lattice results: i. spectral densities from 
MEM analysis of the correlators; ii. static quark free 
energies F{R); iii. quasiparticle masses; iv. bulk ther- 
modynamics p{T). We conclude that everything indeed 
fits well together, provided the role of multiple (colored) 
bound states is recognised and included. 

We have parameterized recent lattice data on free en- 
ergies for static quarks, calculated the corresponding ef- 
fective potentials and solved the Klein-Gordon equation 
for charmonium, light quarks and singlet gg cases. We 
have found that the bound states exists at T below some 
calculated zero-binding points. Our reported range of 
temperatures for charmonium and light mesons above Tc 
agree rather well with what was seen from lattice corre- 
lators using the MEM [15]. 

Our studies of relativistic effects have found that the 
systems in question are not yet very relativistic, so that 
the relativistic corrections to the potentials do not exceed 
10%. They have minor effect on binding, but are more 
important for the wave function at the origin. Like the 
authors in [22] we concluded that some non-perturbative 
interactions for light quarks, completely absent for static 
ones, should exist in order to bring the pion and sigma 
mass to zero at T = Tc. It is believed to be due 
to instanton-antiinstanton molecules, and is quasi-local: 
thus one can treat this interaction like a delta-function 
potential, using the wave function at the origin calculated 
without it. 

Finally, we have assessed the contribution of all these 
bound states to the bulk pressure of the system. We have 
shown that as the bound state inches towards its end- 
point, its contribution to the pressure becomes partially 
compensated by a repulsive effective interaction between 
the unbound quasiparticles. The contribution of the vir- 
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tual level above zero quickly disappears^. All in all, we 
have found that the summed contributions from the large 
number of colored bound/free states above agree well 
with the bulk lattice pressure. 
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APPENDIX: 
THE KLEIN-GORDON EQUATION FOR A 
COULOMB PLUS QUASI-LOCAL POTENTIAL 

In this appendix we will use notations for a single body 
KG equation, in schematic notations (E—V)"^ —p^ = Af^, 
for a particle of mass M in a potential V, chosen to be a 
superposition of a screened Coulomb and an additional 
local term 



V 



a 



(r + 0.001/M) 



exp(-M£,r) - Ud{r) 



(72) 



where the 6{r) is the "nonlocal delta function" . For rea- 
sons related with instantons we will use it in the form 



US{r) = U 



1 



(r2 + p2) 



,2^3 



(73) 



with the size parameter p to be chosen below to be 
p = 1/(0. 6M). Note that we have chosen to ignore the 
coupling constant running, and regulate the Coulomb sin- 
gularity. 

Together with the quasiparticle mass(es) M we 
thus already have four parameters a, Md,U, M . The 
relevant ones are the 3 dimensionless combinations 
a, Mjy/M, U/M'^. To get familiar with this problem, we 
first studied the geometry of the fixed energy surfaces, 
E — const, in the parameter space. A section of those 
surfaces with two coordinate planes is shown in Fig. 7. 

In Fig. (a) one can see that the (regularized)** Coulomb 
and the instanton-induced potential are kind of comple- 
mentary to each other, except near the origin: Coulomb 
always have levels why the quasi-local potential does not 
for U < Umin ~ 20. On another plane, as AId grows and 



^All of this happens gradually like in the ionoization pro- 
cess, so that the disappearence of a state does not lead to 
singularities and phase transitions of any kind. 

** As it is known from the exact solution of the unregularized 
KG equation, the solution gets singular at a = 1/2. There- 
fore all our results for a > 1/2 are actually sensitive to the 
regularization used. 



screen the Coulomb field, one needs stronger coupling to 
keep the same binding. 





I I I I I I I 1 1 I I I 1 1 I I I 1 1 I I I I 1 1 I I I 1 1 I I I I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 




FIG. 9. A schematic plot explaining the geometry of the 
surfaces of fixed binding in the 3-dimensional parameter 
space. Its projection on the a, U /M^ plane and on the 
a^Mo/M plane are shown in two subsequent plots. In 
the former case the values of the energy (in units of the 
mass) is E/M = 0, 0.5, 0.75, 0.9, 0.96, 0.99, from top to bot- 
tom. In the latter case only the lines corresponding to 
E/M = 0.5,0.9,0.99 are shown (right to left). 

The very left line, corresponding to E/AI — 0.99 or 
only 1% binding, is close to the "zero binding line" (ex- 
cept that it reaches the origin a = 0, Md = 0), to the left 
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of which the potential in question has no bound states at 
all. 

Wc have not plotted the third projection, as for a = 
the value of Md is irrelevant and all lines depend on the 
U value only. The zero binding line, separating the un- 
bound states from the bound ones, starts at the value 
of Umin ~ 20 already mentioned. Combining 3 projec- 
tions for the same binding, one can now well imagine the 
location and the shape of all constant energy surfaces. 
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